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Abstract 

A classification, according to invariant theory, of non-constant invariant Abel ODEs known 
as solvable and found in the literature is presented. A set of new integrable classes depending 
on one or no parameters, derived from the analysis of the works by Abel, Liouville and Appell 
[H |[ 0, is also shown. Computer algebra routines were developed to solve ODEs members 
of these classes by solving their related equivalence problem. The resulting library permits a 
systematic solving of Abel type ODEs in the Maple symbolic computing environment. 

PROGRAM SUMMARY 

Title of the software package: Extension to the Maple ODEtools package 
Catalogue number: (supplied by Elsevier) 

Software obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland (see application 
form in this issue) 

Licensing provisions: none 

Operating systems under which the program has been tested: UNIX, Macintosh, Windows (95/98/NT). 

Programming language used: Maple V Release 5 

Memory required to execute with typical data: 16 Megabytes. 

Keywords: Abel type first order ordinary differential equations (ODEs), equivalence problem, integrable 
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Nature of mathematical problem 

Analytical solving of Abel type first order ODEs having non-constant invariant. 
Methods of solution 

Solving the equivalence problem between a given ODE and representatives of a set of non-constant invariant 
Abel ODE classes for which solutions are available. 

Restrictions concerning the complexity of the problem 

The computational routines presented work when the input ODE belongs to one of the Abel classes considered 
in this work. This set of Abel classes can be extended, but there are classes - depending intrinsically on many 
parameters - for which the solution of the equivalence problem, as presented here, may lead to large and 
therefore untractable expressions. When the invariants of a given Abel ODE depend on analytic functions, 
the success of the routines depends on Maple's ability to normalize these invariants and recognize zeros (this 
is well implemented in Maple, but it may nevertheless not work as expected in some cases). Also, when 
the solution for the class parameter depends on other algebraic symbols entering the ODE being solved, the 
routines can determine this dependency only when it has rational form. 

Typical running time 

The methods being presented here have been implemented in the framework of the ODEtools Maple package. 
On the average, over Kamke's [Q first order Abel examples (see sec. the ODE-solver of ODEtools is now 
spending w 6 sec. per ODE when successful^ and « 11 sec. when unsuccessful. The timings in this paper 
were obtained using Maple R5 on a Pentium II 400 - 128 Mb. of RAM - running Windows98. 

Unusual features of the program 

These computational routines are able - in principle - to integrate the infinitely many members of all the non- 
constant invariant Abel ODE classes considered in this work. Concretely, when a given Abel ODE belongs 
to one of these classes, the routines can determine this fact, by solving the related equivalence problem, and 
then use that information to return a closed form solution without requiring further participation from the 
user. The ODE families that are covered include, as particular cases, all the Abel solvable cases presented in 
Kamke's and Murphy's books, as well as the Abel ODEs member of other classes not previously presented 
in the literature to the best of our knowledge. After incorporating the new routines, the ODE solver of the 
ODEtools package succeeds in solving 97 % of Kamke's first order examples. 



LONG WRITE-UP 



1 Introduction 

From some point of view, after Riccati type ODEs, the simplest first order ordinary differential 
equations (ODEs) are those having as right hand side (RHS) a third degree polynomial in the 
dependent variable, also called Abel type ODEs|^ 

y' = f3y^ + f2y' + fiy + fo (i) 

where y = y{x), and /o, /i, /2 and /s are analytic functions of x. As opposed to Riccati ODEs, 
for which integration strategies can be built around their equivalence to second order linear ODEs, 
Abel ODEs admit just a few available integration strategies, most of them based on the pioneering 
works by Abel, Liouville and Appell around 100 years ago ||2|, ^, In those works it was shown 
that Abel ODEs can be organized into equivalence classes. Two Abel ODEs are defined to be 
equivalent if one can be obtained from the other through the transformation 

{x = F{t), y{x) = P{t)u{t)+Q{t)} (2) 

where t and u{t) are respectively the new independent and dependent variables, and F, P and Q 
are arbitrary functions of t satisfying F' ^ and P ^ 0. 

Integration strategies were then discussed in |^, around objects invariant under Eq.(||)0 
(herein called the invariants) which can be built with the coefficients {/s, /2, /i, /o} and their 
derivatives. To each class there corresponds a different set of values of these invariants, and actually 
any one of them (we shall pick one and call it the invariant) is enough to characterize a class. 
A simple integrable case happens when the invariant is constan^ the solution to the ODE then 
follows straightforwardly in terms of quadratures, as explained in textbooks ^. On the contrary, 
when the invariant is not constant, just a few integrable cases are known and the formulation of 
solving strategies based on the equivalence between two such Abel ODEs, one of which is integrable, 
appears to be only partially explored in the literature, and not explored in general in computer 
algebra systems. 

Having this in mind, this paper concerns Abel ODEs with non-constant invariant and presents: 

1. A discussion and classification of the integrable Abel ODEs found both in Kamke's book and 
in the works from the late 19*'* and early 20^^ century by Abel, Appell, Liouville and other 
sources; 

2. A set of new integrable Abel ODE classes - some depending on arbitrary parameters - derived 
from those aforementioned works; 

3. An explicit method of verifying or denying the equivalence between two given Abel ODEs, 
one of which we know how to solve since it represents one of the above mentioned classes; 
and in the positive case, a way to determine the function parameters F, P and Q of the 
transformation Eq.(^) which maps one into the other; 

^For convenience, in this work, by "Abel ODEs" we mean Abel ODEs of first kind, since Abel ODEs of second 
kind can always be transformed into first kind by a simple change of variables. 
^The invariants change in form for F{t) ^ t, but keep their value. See Eq.(|^. 
^There exists one invariant such that if it is constant then the other invariants are as well. 



4. A computational scheme to resolve the equivalence problem in the case of parameterized 
classes, including the determination of the value for the class parameter when the equivalence 
exists; 

5. A set of computer algebra (Maple) routines implementing the algorithms presented in items 
(^) and above, to systematically solve - in principle - any Abel ODE belonging to one 
of the classes, parameterized or not, presented here and for which a closed form solution is 
known (items (||) and (^) above). 

Item (1) is interesting since the Abel ODEs shown in textbooks in general, including Kamke's 
book, are displayed without further classification, and in fact many of them belong to the same 
class. This classification in terms of invariant theory is also necessary in a computational scheme 
for solving Abel ODEs as the one being presented, and we have not found it in other references. 
The integrable classes mentioned in (^ are new to the best of our knowledge, although directly 
or indirectly derived from previous works. The formulation of the equivalence problem mentioned 
in (^) is the one given by Liouville in |^], is systematic and does not involve solving any auxiliary 
differential equations^. Concerning item @), the idea can be viewed as a way of avoiding the 
untractable expressions which one would encounter when making direct use of Liouville's strategy 
with parameterized classes. The strategy presented is applicable when there exists a solution for 
some numerical values of the parameter, or when this parameter is a rational function of other 
symbols entering the input ODE. Regarding item (^), the implementation presented here is, as 
far as we know, unique in computer algebra systems in its ability to solve non-constant invariant, 
parameterized or not, Abel ODE classes. 

The paper is organized as follows. In sec. |2|, the basic definitions and the classic formulation of 
the equivalence problem in terms of invariants is reviewed and shown to apply straightforwardly to 
the case of a non-parameterized class. In sec. ^ it is shown how these ideas can be complemented 
by taking advantage of computers to tackle the equivalence problem in the case of a parameterized 
class. Section ^ presents a classification of the integrable classes we have found in the literature 
with some additional comments as to their derivation. In sec. ^ new integrable Abel classes are 
presented. In sec. |^, a test-suite for the routines presented is discussed and statistics are shown 
describing the performance with this test-suite as well as with Kamke's first order examples. Finally, 
the conclusions contain some general remarks about this work and its possible extensions. 

Additionally, we present in the Appendix a table listing the distinct Abel ODE classes that we 
have found, representative ODEs from each class, and their respective solutions. 



2 Classical Theory for Abel ODEs 

In general, Abel type ODEs can be studied using two related concepts: invariants and ODE 
equivalence classes. We define two Abel ODEs to be equivalent^ if one can be obtained from the 
other using a transformation of the form Eq.(§). The equivalence class containing a given ODE 
is then the set of all the ODEs equivalent to the given one. We note that although the infinitely 
many members of a class can be mapped into each other by using Eq.(|2|), there are also infinitely 
many disjoint Abel classes (Eq.(^) is not sufficient to map any Abel ODE into a given one). 

To each class one can associate an infinite sequence of absolute invariants To see this, con- 

sider two Abel ODEs, the first Eq.(|l[), the second obtained from Eq.(||) through the transformation 

*An approach somewhat similar to this one by Liouville is discussed in [^. 
^For a more formal definition of class see |7f 



Eq.(l) 

u' = f3u'' + f2U^ + fiu + fo (3) 
where the coefficients {/o, /i, /2, /a}, are related to the those of Eq.(|I]) by 



/o 



F' {MF) + h{F) Q + h{F) Q2 + /3(F) Q3) _ Q' 



P 

P' 



fi = -p-F'(h{F) + 2f2{F)Q + 3MF)Q' 
h = PF'(/2(F)+3/3(F)q) 

/a = P^F'hiF) (4) 

Following 0, we call an absolute invariant oi'EiC{.@) a function /(/, x) of the coefficients {/o, /i, /2, /s} 
and their derivatives with respect to x such that, for all {F, P, Q} in Eq.(^), 

Hf,t)\f=f(f,t) = Hf,x)\x=F{t) (5) 

where / = f{f,t) represents the coefficients {/o, /i, /2, /s} and their derivatives with respect to t, 
expressed using Eq.(Q). 

Similarly, we call a relative invariant a function, say s, of the coefficients of Eq.(|l]) and their 
derivatives such that when changing variables using Eq.(^, the resulting expression is equal to the 
original one up to a factor, say ips, dependent uniquely on the functions F, P, and Q in Eq.(^) and 
independent of the coefficients themselves Q: 

■<f)\f=fif,t) = MF, P, Q) s(/)U=F(t) (6) 
Liouville showed that in the case of Abel equations there is a relative invariant of weight 3 

S3 ^ Ml + \ - /1/2/3 + hf2 - f2f3^ (7) 

which can be used to recursively generate an infinite sequence of relative invariants S2m+i of weights 
2m + 1 respectively^, through the formula 

S2m+1 = h S2m-1 " (2 m - 1) S2m-l {^f^ + /l /s " (8) 

As is clear from this definition, the product of two relative invariants respectively of weights n and 
m is a relative invariant of weight n + m, and by dividing any two relative invariants of equal weight 
one can generate an infinite sequence of absolute invariants 

h = ^, h = ^3 = — , etc... (9) 



^In the case of S3, 1^33 = (F'P)^; the weight n refers to the degree of <^s„ with respect to (F'P). 



In [Q], Appell showed that this sequence can also be obtained from two basic absolute invariants 
- say Ji, Jo, by expressing Ji as a function of Jq and then differentiating the result with respect 
to Jq. As a consequence, if Ii is constant then all the other ones will be too. This fact allows one 
to identify the constant character of the invariants in Eq.@ by looking at just the first one. We 
note also that there are infinitely many different classes having Ii constant, related to the infinitely 
many possible constant values Ii can have. 



2.1 Integration strategy 

A description of a method of integration when the invariants are constant]^ is found in the works by 
Abel Q, Liouville and Appell |^. In such a case, all members of the class can be systematically 
mapped into a separable first order ODE (representative of the class), by appropriately choosing 
F, P and Q in Eq.(P); see for instance and 

A quite different situation happens when Ii is not constant. In such a case, relatively few 
integrable Abel ODEs are known, and the integration methods used to solve each of them depend 
in an essential way on non- invariant properties of the coefficients /. Those methods are then useless 
for solving the other infinitely many members of the same classes, unless one can solve the related 
equivalence problems; i.e., determining - when they exist -the values of F, P and Q in Eq.(^) linking 
two Abel ODE which belong to the same class. 



2.2 Identifying an ODE as member of a given Abel ODE class 

Consider two Abel ODEs; the first one given by Eq.(|l|), and a second one being of the same form, 
but with coefficients /o, /i, /2 and /a. The problem now is to determine whether the second 
Abel ODE can be obtained from Eq.(||) by changing variables using Eq.(|2|). This problem can be 
formulated by equating the coefficients between the transformed equation, obtained by applying the 
transformation Eq.@ to Eq.(||), and the second Abel ODE, resulting in Eq.(^), which can be seen 
as an ODE system for {F, P, Q}. To solve this system, following Liouville |p, we first note that 
the absolute invariants corresponding to the two Abel ODEs don't depend on P or Q (see previous 
section). Hence the function F entering Eq.(^) can be obtained by just running an elimination 
process using two of these absolute invariants, for instance Ii = sf/sl and I2 = s^sy/s^: 

"-75-3!^=^*) ^ - ^ - —4-\x=F{t) (10) 
S3 63 S3 63 

As discussed in ^, the existence of a common solution F{t) to both equations above (such that 
F' 7^ 0) is the necessary and sufficient condition for the existence of a transformation Eq.(^) relating 
the two Abel ODEs. Once F is known, the system Eq.(Q) becomes trivial in that P and Q can be 
re-expressed in terms of F by performing fairly simple calculations. In the case of interest of this 
work - non-constant invariant^ - the resulting expressions are: 

~ 2 ~ ~ 

D^.x F' fs S3 F' f2 fsS^ - f2f3S3 

^(t) = ,2- \x=F{t) Q[t) = n f 2~ \x=Fit) (11) 

fs S3 3 /s S3 

where {fi, fi} with z : — > 3 are the coefficients of the two Abel equations, S3 is the relative 
invariant Eq.(^) expressed in terms of /, and s~3 = S3|j.^j. ■ 

^In 1^] it is also shown that in the constant invariant case the problem can also be formulated in terms of the 
symmetries of these ODEs. 

*When the invariant is non-constant, S3 7^ 0. 



Concerning the explicit solution F{t) for Eq.(^), we note that our interest in solving the 
equivalence problem is in that it leads directly to the solution of other members of an Abel class, 
when the solution to a representative of the class is known. In turn, all the solvable classes we 
are aware of have a representative with rational coefficients (see the Appendix), and hence also 
rational invariant^. Hence, assuming that one of the two Abel ODEs has rational coefficients and 
that Eq.(|lO|) was obtained by applying Eq.(§) to it, the system Eq. (|lO|) will always be polynomial 
in F{t). In such a case, when a common solution F{t) to both equations exists, the resultant 
between these polynomials will be zero Q; i.e.: there will be a common factor, depending on F and 
t and representing the common solution, which can be obtained by calculating the greatest common 
divisor (GCD) between the two equations in Eq.(|l^). Conversely, if that CCD does not depend on 
F, a transformation Eq.(§) linking the input equation to Eq.(||) does not exist. That the dependence 
on F of this GCD is a necessary condition for the existence of the desired transformation Eq. (|^) is 
a consequence of the validity of Eq. @ and hence the system (|To|) . A proof of its sufficiency was 
given by Appell in Q. 

The whole process just described to determine the equivalence between two given Abel ODEs, 
one of which is rational in x, can be summarized as follows: 



1. Calculate two absolute invariants, set up the system Eq.(pl)[), and calculate the GCD between 
the two equations; 

2. When this GCD does not depend on F, the ODEs don't belong to the same class; otherwise 
determine an explicit expression for F{t) from the result of the GCD calculation; 

3. Plug this value for F into the formulas Eq.(O) to determine the values of P{t) and Q{t), 
arriving in this way at the transformation Eq.(^) mapping one Abel ODE into the other. 

Example: 

Consider the two non-constant invariant Abel ODEs 

,_ {fx- f) (, , , y 



y 



f{x-f)y^ + y^)-y (13) 

V / .1? 



2(/ + 3x) V' ' / X 

where in the above / = f{x) is an analytic (arbitrary) function. As in the typical situation one of 
these ODEs has invariants rational in x and we know its solution; i.e. for Eq.(^) we have 



C, + ^y'^^-^'y-^ + 2 arctan ( , ^ + '^ ) = (14) 
y \y/y^x - Ay - I J 

where Ci is an arbitrary constant. We would then like to determine whether there are functions 
{F, P, Q} so that Eq.(|l^) transforms under Eq.(^) into Eq.(|T3|), and if so, determine {F, P, Q} 
and use them together with Eq.(p^) to build the answer to Eq.(p^). For this purpose, we start 
(step dl])) by computing the relative invariants S3, S5, sy, leading to Eq.(|lO|) 



®0n the other hand, there is no reason to expect that the second of the two Abel ODEs being tested for equivalence 
also has rational invariants. If however the invariants of both Abel ODEs are rational in t and the coefficients are 
numbers, then it is also possible to determine F{t) by performing a rational function decomposition as mentioned in 



(87/4+9/2+184*2) 
t(9/+31i)'^ 



(280+105 F+9F2)'^ 

(81 /3+1431 /2t+7185 /t2+10903 t^){9 /+31 1) 
(87/t+9/2+184t2)2 



(1674 F2+10290 F+81 F3+19600) (40+9 F) 



(15) 



(280+105 F+9F2)2 

where in the above F = F{t) is the function we are looking for and / is taken at x = t. We recall 
that when Eq.(|T^) has non-constant invariant the denominators in above are not zero since 53 7^ 0. 
Calculating the GCD between the numerators of the expressions above (step (|^)) and equating this 
GCD to zero, we obtain^ 



27{t + tF-f) = 
from where the common solution F(t) to both equations is given by 

fit) 



1 



(16) 



(17) 



Substituting this value of F into Eq.(|ri|), a transformation of the form Eq.(^) mapping Eq.(|T^) 
into Eq.([l3|) is finally given by 



{x 



fit) 



1, yix) = tuit)} 



(18) 



from where by changing variables in Eq.(|T^) using the transformation above and renaming the 
variables (t — > x, u ^ y), the solution to Eq.(|l3|) is obtained 



Ci + 



x'^y^ 



Axy — 1 



( 



2 arctan 



xy 



(l + 2xy) 



y (I - 1) a:2y2 _ 4 - 1 ^ 



(19) 



3 Parameterized Abel ODE classes 



We formulate here the equivalence problem in the case of parameterized classes. By "parameterized 
class" we mean an (Abel) ODE class depending on symbolic parameters which cannot be removed 
by changing variables using Eq.(^). The interest in parameterized solvable classes is clear: to each 
set of values of the parameters corresponds - roughly speaking - a different Abel class|^. Hence, 
a formulation of the equivalence problem for parameterized classes enables one to solve all the 
members of infinitely many classes at once. 

In order to simplify the discussion, we consider the problem of an Abel ODE class depending on 
just one parameteiPl, say C. Also, we distinguish between two different types of problems: one is 
when the equivalence problem has a solution for a specific numerical value of C; the other happens 
when to have a solution it is required that C assumes symbolic values, for instance in terms of other 
symbols entering the input ODE. We discuss first the numerical case, and in the next subsection we 
show how the symbolic case can be mapped into many numerical problems - when the parameter 
depends on other symbols in a rational manner - by using rational interpolation methods. 

^"We note that when f{t) is an algebraic mapping involving varied analytic functions, Maple's procedures to 
simplify and put them in normal form may fail and consequently the GCD computation may not be successful. 

There may be particular different sets of values for which the resulting ODEs will nevertheless belong to the same 
class. 

^2The integrable classes presented in the literature depend at most on one parameter (see sec. 0). 



3.1 Solution for some numerical value of C 



To facilitate the exposition we present the discussion around a concrete example. Consider the 
equivalence problem between a given Abel ODE, for instance, 

y' = 8 ^ = ^ + 4 ^ + ^ 20) 

x' X^ X 

and the one presented in Abel's memoires |2| 

(C.' + x^ + l),3^ (21) 

x-^ 

If this equivalence exists, then it exists for a specific value of the parameter C since there is no 
solution for arbitrary C (the existence of such a solution would mean the class does not really 
depend on any parameter). Hence, the common solution F{t) to the system Eq.(^) will not show 
up until the correct value of C is determined^, invalidating the itemized algorithm of the previous 
section. 

A natural alternative to this problem would be to take one more absolute invariant, for instance, 
S3S7/S5, so that our system Eq.(p^) becomes 

n _ «~5^ si S3S7 S3S7 _ S5S7 S5S7 

^ - j5 - ZE\^=F{t) ^-T2 ;r\^=F{t) ^--^4 ;r\^=F{t) [^^) 

S3 63 S5 6g S3 S3 

and search for a solution to this problem such that F' ^ 0, C = 0. For that purpose, eliminate 
C from the first and second expressions above by taking the resultant with respect to C, obtaining 
- say - Ri. In the same way, eliminate C from the first and the third expressions of Eq.(p^) 
obtaining i?2- Hence, when a solution exists, the resultant between Ri and R2 with respect to F 
will vanish. In other words, the algorithm of the previous section will work if instead of performing 
the calculations over the expressions Eq.([To|) we perform them over Ri and R2. The GCD between 
i?i and i?2 will then return the factor depending on both F and t, whose solution is the function 
F{t) we are interested in. 

This method, simple and correct in theory, unfortunately does not work in practice because the 
expressions tend to grow in size so much that the computation of the first of these three resultants 
may not be possible, even with a simple example such as the one shown above. The problem resides 
in the fact that multivariate CCDs and resultants are quite expensive operations for the current 
symbolic computation environments. 

An alternative to this problem consists of reducing it to a sequence of bivariate GCD and 
resultant calculations, for which the available algorithms are relatively fast. The idea can be 
summarized as follows. 

1. From the previous considerations, when a solution to the equivalence problem exists, the 
resultant between any two of the expressions in Eq.(^2|) will not vanish for any value of t, 
since we haven't introduced the correct (unknown at this point) value of C. Hence, if we 
insert in Eq.(p2[) a numerical value^ for t and calculate the GCD between any two of the 
resulting expressions, this GCD cannot contain any factor depending on F. This gives us a 
first "existence condition" test for the solution before proceeding further; 



^''There may be more than one solution C. 



^*We note there may exist "invahd evaluation points"; roughly speaking to avoid thisproblem this evaluation point 
must not cancel any of the coefficients of the variables remaining in the system - see 



2. When Eq.(p^) evaluated at t = number passed the test of the previous step, take two of the 
resulting three expressions and calculate their resultant with respect to F, obtaining, say, Ri. 
Then take a different pair and calculate their resultant with respect to F again, obtaining, 
say, i?2- Neither of these resultants will vanish since the GCD calculations of the previous 
step showed no factor depending on F. Also, the calculation of Ri and i?2 is now quite 
simpler since the expressions do not involve t; 

3. Then if a solution to the problem exists, the GCD between Ri and R2 will yield a factor 
depending on C; equating it to zero and solving it for C will give the common solution C for 
Ri and i?2- More precisely, what we will get in this way is a set of candidates (including 
among them the correct value) for C; not all of them will necessarily lead to a solution F{t) 
to the original problem; 

4. We now plug these candidates for C into Eq.(^), one at a time, receiving a system of three 
expressions involving again only two unknowns, now F and t. If there is a common solution 
F(t) to these expressions, the resultant with respect to F between any two of them will vanish. 
Hence, the GCD between those two expressions will contain a factor depending both on F 
and t; equating this factor to zero and solving for F leads to the solution F{t). 

Returning to our example of determining the equivalence between Eq.(|20|) and Eq.(^), the itemized 
procedure just outlined runs as follows. 

According to step (||), t = is tried first, but it is found to be an invalid evaluation point. The 
next value of t to try, t = 1 turns out to be valid, so Eq.(p2|) was evaluated at i = 1; the GCDs 
between any two of the three resulting expressions do not depend on F, so this first test for the 
"existence" of a solution passed. 

Continuing with step (|2|), the calculation of Ri and R2 is performed without problems concerning 
the size of the expressions. 

The GCD of step (^) results in the three factors: 36C — 5, C + 1 and 9C — 2; equating them to 
zero and solving them for C we arrive at three candidates for C. 

In step (^), plugging each of these candidates one at a time into Eq.(^) and taking the GCD 
between two of the three resulting expressions we note that C = 5/36 does not lead to any factor 
depending on both F and t, but C = — 1 leads to such factor: F^f^ — 1. So that for C = — 1 the 
problem admits two solutions: F = itl/t^. 



Finally, by introducing F = into the formulas for P and Q Eq.(ll) we arrive at the 
transformation of the form Eq.(|2|) mapping Eq.(|2lD into Eq.([20|) 

u{t) 



^,y = -2^\ (23) 



and hence by applying the same change of variables to the answer of Eq.(]2l|) and substituting 
C = — 1 we obtain the answer to Eq. (pO|) . 

A remark however is in order: if, in step (^ of the algorithm just described, the numeric 
candidates for C involve fractional powers of rational numbers, the Maple system may then enter not 
efficient expensive computations, exhausting the system resources before determining the expression 
F{t) solving the problem in step (Q). The root of this limitation seems to be in the absence in 
Maple of built-in normalization for such "numeric radicals"]^. 

built-in normalization of radicals is implemented in the computer algebra system "Mathematica" . For typical 
problems we tried where Maple exhausted the system resources trying to determine the solution F{t), we exported 



3.2 Solution when the parameter C is some rational function of other symbols 



When C assumes symbolic values, for instance it depends on other symbols - say {a} - entering the 
input ODE, if this dependency is rational it is possible to map the determination of C{a) into a 
sequence of problems having for solution a numerical value of C. In turn each of these numerical 
problems can be tackled using the algorithm of the previous subsection. The idea consists of 
attributing numerical values to the symbols {a} entering the invariants in order to determine C{a) 
by means of a rational interpolation. To simplify the presentation we first discuss the case when 
{a} consists of a single parameter, and then show how to extend the algorithm to the case in which 
{a} consists of many parameters by commenting on a concrete example. So when a consists of just 
one parameter this rational interpolation scheme is summarized as follows: 



1. Take the system Eq.(22) and attribute a numerical value to a (check for possible wrong 



evaluation points), so the resulting system depends on just x, C and F; 

2. Enter step (|l]) in the algorithm of the previous subsection and run all the steps: 

(a) If there is no solution for C and F such that C = and F' ^ then quit the process 
- the input ODE does not belong to this class, or the solution involves a non rational 
dependency of C on {a}. 

(b) If however a solution for C and F was found, record the values for C and a - they 
represent a point of a curve C(a); 

3. Using the points recorded so far, interpolate C as a function of a and test if this interpolated 
value already solves the problem for an arbitrary of^: 

(a) If so, the problem has been solved; 

(b) Otherwise change the evaluation point of a and re-enter step (|2|) of this enumeration; 

We note that the rational interpolation of C(a) requires the knowledge a priori of the polynomial 
degrees in a of both numerator and denominator. That information is not available in advance, but 
we know that these two degrees sum to n — 1, where n is the number of points being interpolated. 
So, when performing the test in step (|3|) and before going to step ( pb| ) we actually test all possible 
different interpolations, starting with the maximum possible degree for the numerator and finishing 
with the maximum possible degree for the denominator. 

Concerning the extension of this algorithm for the case when {a} involves more than one pa- 
rameter, for instance {a, b}, this extension is easy and better illustrated with an example. Consider 
the equivalence problem between the Abel ODEs 



, C (2x^C -2)y^ -SCy"^ + Cxy , , a (2x'^a - 2 b) y^ - Say'^b + abxy 
1 — x^C b{b — x^a) 

which solution is just the identity x ^ x, y ^ y , but only exists when C = | (we choose the 
identity without loss of generality and so that the solving process is easy to follow). We want to 
determine C = | (rational function of two parameters) using the algorithm just described. 

the mathematical expressions involving radicals to Mathematica and noticed that F{t) was determined in this other 
computer algebra system in reasonable time. 

^''This is done by restoring the symbol "a" in the system obtained in step (□) and checking if the system is satisfied. 



We start with step (|l|), take the system Eq.(^) corresponding to Eq.(^) and attribute numerical 
values for the first parameter, a, checking for possible wrong evaluation points - we end up evaluating 
this system at a = 1. However, the system still depends on the second parameter, 6, so we attribute 
numerical values to b too - check for possible wrong evaluation points - and hence end up evaluating 
the system altogether at a = 1, 6 = 2. The resulting system now only depends on x, F and C. 

So we enter step (§) (actually the whole algorithm of sec. with this evaluated system of 
the previous step and detect that C = 1/2 leads to the solution F = x; so a solution C' = and 
F' / exists. 

Hence, in step (|3|) we interpolate C{b) (so far we have just one point, so C{b) = 1/2) and test 
this value for arbitrary b, verifying that the interpolation is still incomplete: the system Eq.(|2^) 
with a = 1 and arbitrary b is not satisfied by C = 1/2, F = x. 

We are then in step (|^.^); so attribute a new numerical value to 6, hence evaluate the system 
Eq.(22) at a = 1, 6 = 3 and re-enter step (§) finding that C = 1/3 leads to F = x. 

We record this new point (step (3.|2b|)) and interpolate - now using the two points obtained so 
far - our first interpolation is of degree 1 in b: C{b) = (5 — b)/6. This interpolation however does 
not solve the problem for arbitrary b. So we increase by one the degree in b of the denominator in 
the interpolation, resulting in C{b) = 1/6, and verify that this second interpolation indeed solves 
the system for arbitrary 6 - so the interpolation for 6 at a = 1 is complete. 

We return then to the interpolation of C{a, b) with respect to a, record that for a = 1 there is 
a solution C{b) = 1/6, and test this solution for arbitrary a (step (P) with respect to a), verifying 
that the interpolation for a is still incomplete. 

Hence we re-enter step (||) with new evaluation points a = 2, 6 = 4, and re-start the process of 
determining C(6). 

Following the same steps just described, for a = 2 we find that C(6) = 2/6 solves the problem 
for arbitrary 6, leading to the second point in the interpolation of C{a, 6) with respect to a. 

So we are now in step (^) again, and interpolating C using the two points obtained so far we 
find C(a, 6) = a/6 - this value of C(a, 6) is verified to satisfy the system Eq.(^) for arbitrary a, and 
thus the problem has been solved. 

The algorithm just described, though expensive in computations, is successful in solving the 
equivalence problem in reasonable time for typical situations (see sec. |). In this example, for 
instance, it took 12 seconds to: setup the invariants and the system Eq.(22|), simplify this system 
to a normal form, run all the items of the algorithm just described to determine C{a, 6) and F{t), 
then determine P{t) and Q{t) according to Eq.([rT|), and finally return a solution to the second of 
the ODEs in Eq.(p|). 



3.2.1 Remark on the existence of multiple solutions for the class parameter C(q) 

The interpolation algorithm just described is valid provided there is only one solution curve C(a); 
otherwise we may end up trying to interpolate C{a) using points which belong to different solution 
curves, leading nowhere. In turn, the existence of many curves C{a) solving a given problem is 
related to the existence of symmetries in the invariants (of the ODE representative of the class 
we want to match) entering Eg . (p2|) . Concretely, if the mapping {C — > k{C), x — > 4>{x,C)} is a 
symmetry of these invariants, then if {C, F{x)} leads to a solution for the equivalence problem, 
consequently {k{C), (j){F[x),C) } will also lead to a (different) solution. A concrete example of this 
situation is discussed in sec. |6.1| . 

Concerning detecting this situation, we note that if the mapping {C k{C), x (j){x,C)} is a 
symmetry of the invariants, then the inverse mapping is also a symmetry, and since these invariants 



are rational in both x and C, the form of such a symmetry mapping is 



aC + b f{C)x + g{C) 

that is, a fractional hnear mapping, where {a, b, c, d} and {/(C), 5(C), /i(C), j(C)} are in principle 
constants and functions to be determined. Regarding the use of interpolation methods, the problem 
happens when k{C) ^ C. A first manner of detecting this problem then consists of setting up the 
system of equations and inequations 



/(x,C) =/(x,C) 'P=r^ F,, = 0, C = 0, Ac(C)/C (26) 

where I{x,C) = sf/sl is the first invariant and F, G and k are the unknowns to be determined, 
and seeing if this system is consistent. This check for consistency can be performed by simplifying 
this system with respect to its integrability conditions - for this purpose we used the diffalg [O] 



and RIF [14| Maple packages. When the system has no solution, this fact is detected by these 
packages; otherwise the related symmetry of the invariants is obtained from the output of these 
packages directly. 



4 Integrable Abel ODE classes found in the literature 

This section is devoted to a compilation of integrable Abel ODE classes found in the literature. 
The compilation is not intended to be complete, but it nevertheless covers various of the usual 
references; mainly Kamke's and Murphy's books [||, |5|, and the original works by Abel, Liouville 
and others on these subjects |^, ^ ^ [l0| ]. 

One of the noticeable things in these references is that the presentation of integrable cases lacks 
a classification in terms of their invariants. Consequently, many of these ODEs can actually be 
obtained from one another by means of Eq.(|^, that is, they belong to the same class. Since part of 
this work consisted in writing computer routines addressing the equivalence problem, we performed 
this classification, and therefore present here a more compact collection of integrable Abel ODE 
classes, as opposed to just integrable ODEs. Classes not depending on parameters are labelled by 
numbers (e.g.. Class 1), while those depending on parameters are labelled with letters (e.g.. Class 
A) . The solutions to the representatives of these classes are presented altogether in a table in the 
Appendix. 

While revising the related literature we also noticed that various of the cases presented in 
books or papers are in fact particular cases of the integrable classes presented by Abel, Liouville 
and Appell in [^, ^, |^. In turn the methods they used to obtain new integrable classes seem to 
be forgotten or not mentioned elsewhere. So, it appeared reasonable to start by reviewing and 
analyzing selected parts of those works in this section, and then show in the next section how, 
starting from these ideas, additional integrable classes can be obtained. 

The first large presentation of integrable cases is due to Abel himself in H]. His idea was to 
consider integrating factors of the form 

^^ = e^(^'^) (27) 
for "Abel" equations written in terms of two arbitrary functions p and q as: 



$ = yy' + p{x) + q'{x)y = 



(28) 



The first non-trivial case discussed in |^ was found by taking r(x,y) as quadratic in y: 



where a, (3 and 7 are arbitrary functions of x. Abel formulated this problem by applying Euler's 
operator to the total derivative obtaining a system easily solvable for a, [3, 7 and p. The 
resulting Abel family has non-constant invariant and is shown in Abel's memoires as depending on 
one arbitrary function q{x) and two arbitrary constants Cj: 

(for the corresponding integrating factor see 0]). Now, for the purpose of building computer 
routines addressing the equivalence problem, it is crucial to determine whether or not a given class 
depends on parameters since, as explained in sec. ^, in such a case the formulation of that problem 
is much more difficult. In the case of Eq.(^), the two parameters Q and the function q{x), can 
be removed by first converting the ODE to first kind using y{x) = l/f(x), and then employing a 
transformation of the form Eq.(^): {x = F{t),v{x) = u{t)^/—2Cl}, with F implicitly defined by 
2Ciq{F) — ty/—2Ci + C2 = 0, arriving at a representative of the class simpler than Eq.(pg|), 

y' = - + y' (30) 

X 



and showing that this class does not depend on parameters. It is then easy to verify that Eq.(|3[ 
is a particular case of a parameterized clast^ derived from Appell's work Q. 

The next integrable case shown by Abel is obtained by considering for Eq.(28) an integrating 



factor of the form fj, = exp (l/(a + /3y)). Proceeding as in the previous case, Abel arrived at 
another integrable ODE class with non-constant invariant, which however (see Q) is a particular 
member of the parameterized class Eq.(^) shown by Abel in the same paper. 
Constant Invariant case 

Abel then considered an integrating factor of the form /i = (a + /3 y)". This ansatz does not lead to 
a non-constant invariant family. However, this is the first presentation we have found of a method 
for the constant invariant case. Liouville, and others after him, rediscovered this method, presented 
in Kamke as due to M. Chini Q, and in Murphy's book as a change of variables mapping an Abel 
ODE into a separable one. A recent discussion of the symmetries of this constant invariant problem 
is found in |^. 

Class "A " depending on one arbitrary parameter 
The next ansatz considered by Abel was 

^^ = {A + yr{B + y)'y (31) 

where A{x) and B{x) are arbitrary functions and a and b are arbitrary constants. By taking b = —a 
Abel showed that a tractable integrable case results: 



yy' + ^((g + 2^) -!i^\+q'y = (32) 



The arbitrary function q{x) can be removed together with the constant Ci by rewriting this ODE 
in first kind format, and then appropriately choosing {F, P, Q} in Eq.(pP; so that a simpler repre- 
sentative of this class depending on only one parameter "a" , is given bj- 



Eq.(|30[) is obtained from Eg. (p8[) taking C = and changing variables {x — it,y — iu{t)}. 

representative of the same class of Eq.js^) is shown in as y' = (^{^^ ~ '^^^^ + ^3" 



y'= + ^ + + (33) 

Class 1 

In ll^, Halphen noted a connection between doubly-periodic elliptic functions and the Abel 
type ODE 

x(8y-l) ^ ^ 

which transforms into itself under infinitely many rational changes of variables, from where he was 
able to determine both a parametric and an algebraic solution for it (see the Appendix). 
Class 2 

In a paper by Liouville mostly dedicated to Abel equations, he discussed the integrable cases 
known at that time (1903), and presented some new ones. Liouville reviewed Abel's work and 
considered for Eq.(p8[) an integrating factor of the form Eq.(^) with r{x^ y) cubic in y, arriving at 
the integrable family y' = 6a a y^, depending on a parameter a. This parameter however can 

be removed by changing variables as in arriving at the integrable 

class free of parameters represented by 

y' = -2y^x + y^ (35) 
Class "B" depending on one arbitrary parameter 



As a generalization of Eq.(|35|), in y] Liouville also presented the parameterized family 

y' + (Smx"^ + 4m'^x + n)y^ + 3xy'^ = (36) 



written in terms of two parameters m and n and which can be mapped into a Riccati ODE 
solvable in terms of special functions. Eq.(^5|) is a member of the class represented by Eq.(|3^) 
after setting m = 0. However, when m = 0, n can be removed from Eq.(|3^) by changing 
variables ^x = t^/n, y{x) = t n(t)/n^/^|, leading to a class without parameters - actually rep- 
resented by Eq. (|35D . In turn, when m 7^ 0, m and n can be "merged" by changing variables 
{y = u{t)/m?',x = mt] and introducing a new parameter a = n/m?, resulting in 



V' 



-(3x2 + 4x + a)y^-32;y2 (37) 
In summary, Eq.(|3^) is not a full 2-parameter class, but instead two classes represented by Eqs.dS 



and (37), respectively depending on zero and one parameters. A simpler representative for this class 
and its solution are found in the Appendix. 

Class 3 

Still in [^J Liouville pointed out that by interchanging the role between the dependent and 
independent variables in Eq.(^5|) one arrives at a different Abel integrable class. After rewriting 
this resulting ODE in first kind format and performing a change of variables of the form Eq.(§), a 
simpler representative of this integrable class is given by 



4.1 Integrable Abel ODE classes shown in Kamke and some others books 



One of the most well known collection of (69) Abel ODEs is the one shown in Kamke's book. This 
collection however makes no distinction between constant or non-constant invariant cases, presents 
ODEs of the same class as different, and does not discuss what would be the representative for 
each class depending on the least number of parameters. A first classification for these Abel ODEs 
is then given by|^: 



Classification 


ODE numbers as in Kamke's book 


4 are too general 


50, 219, 250, 269 


40 constant invariant 


38, 41, 46, 49, 51, 188, 204, 213, 214, 215, 216, 218, 221, 222, 223, 224, 225, 
226, 227, 228, 229, 231, 236, 238, 239, 243, 244, 245, 246, 247, 248, 249, 251, 
252, 254, 255, 260, 261, 262, 264 


24 non-constant invariant 


36, 37, 40, 42, 43, 45, 47, 48, 111, 145, 146, 147, 151, 169, 185, 203, 205, 206, 
234, 235, 237, 253, 257, 265 


10 shown without solution 


40, 47, 48, 203, 205, 206, 234, 237, 253, 265 



Table 1. First classification for the 69 Abel ODEs shown in Kamke's book. 



As mentioned, all constant invariant ODEs can systematically be transformed into separable ODEs 
(see for instance Murphy's book), so that the interesting subset is the one comprising 24 ODEs 
having non-constant invariants. We note also that 10 of these 24 ODEs are shown in the book 
without a solution, and in fact we were unable to solve any of 203, 205, 206, 234, 253 or 265, so 
that the number of integrable cases for us is 18. 

Prom these 18 ODEs (and hence from the 69 Abel type Kamke's examples), only four - those 
numbered: 47, 185, 235 and 237- would really lead to additional integrable classes with respect to 
those presented in the works by Abel, Liouville and Appell. We note however that the examples 47, 
185 and 237 are all members of Class "C" (see Eq.(^)), which can be derived from the work by Abel 

- even when it was not presented in the original work. So that the number of additional integrable 
classes presented in Kamke reduces to one, represented by the example 235. The classification and 
details are as follows. 
Class 4 



{xy + a)y' + by = (39) 

This ODE (K 1.235) is presented in Kamke in terms of two arbitrary parameters {a, 6}; then, a 

change of variables which transforms it into a linear ODE is shown. A simpler representative of 

this class - not depending on parameters - can be obtained by rewriting this equation in first kind 

1 

tvJTj 



format via {x = t,y = jX^ — |} and then changing variables {x = j^,y = ^^^}, leading to 



y' = y'- ^^y^ (40) 

X 

Comments on Kamke's example ^7 

y'-a{x^-x)y'-y' = ^ (41) 



For the ODE 

.2 



^^In this classification, by "too general" we mean: these ODEs cannot be solved without restricting the example 
to a concrete particular case. We excluded ODEs - like those numbered 230 and 232 - which are already of Bernoulli 
type. We note also that the ODEs shown in Kamke without solution can all be transformed into an Emden type 
second order ODE shown in Kamke as 6.74, for which only a general discussion is presented. In turn, a detailed 
discussion on the integrable cases of Emden type ODEs is found in . 



presented in Kamke as K 1.47, there is no solution shown in the book, but instead a suggestion 
of transforming the ODE into a second order one. We followed that suggestion and then ran 
a symmetry analysis, noticing that the resulting ODE will have two point symmetries if either 
= ~ ()2f!^n+-n'^ } or {n = 2, a = ^}, leading to two integrable classes not shown in the book. In 



the former case, from Eq.(^), we arrive at 

, (2n + 2)(x" 



x) 



+ + 6n 



However, this ODE can be transformed into Eq.(^9|) by changing variables {x = 1 1-" , y = 
followed by n = so that it belongs to Class C. In the same line, taking {n = 2, 



(42) 



Eq.(pT|), and changing variables {q 



y = 5/2u{t)t } one arrives at Eq.([49D with a 



6, so 



that this second branch of Eq.(^) is also a member of Class C. 
Comments on Kamke's example 237 



x{y + a)y' + by + cx = 



(43) 



This ODE (K 1.237) depending on three arbitrary parameters {a, b, c}, is presented in the book 
without a solution. We note however that changing {x ^ y, y — > x} leads to an ODE also of 
Abel type and in second kind format. Converting the latter to first kind format via {x = t,y = 
^J^^^ — replacing y ^ y' and running a symmetry analysis, the resulting second order ODE 
has two symmetries when a = —2b, leading to an integrable case. Introducing a = —2b into 



Eq.(^3|), rewriting it in first kind format via {x = t,y 



tu(t) 



+ 2 6} and changing variables 



{x = —^-^7^,2/ = b^t+I) ^ ^^^^^ to ^ simpler representative of the class not depending on any 
parameters: 



2 (x + 4) 



(44) 



However, by changing variables {x = 4 (1 — t ) /t , y = —u{t)t/2} one arrives at Eg. ([4 91) again, this 
time with a = —1/2, so that Eq.(^) is also member of Class C. 

A classification for all these 18 non-constant invariant Kamke examples is then as follows^ 



Class 2 


Class 3 


Class 4 


Class A 


Class B 


Class C 


Class D 


36, 40 


145, 147 


235 


257 


42, 43 


45, 47, 48, 151, 185, 237 


37, 111, 146, 169 



Table 2. Classification for the 18 non-constant invariant solvable Abel ODEs in Kamke's book. 

where classes C and D are defined in sec. |. In summary, all but one of Kamke's 58 solvable 
examples (18 non-constant invariant + 40 constant invariant) are particular cases of the integrable 
classes presented by Abel, Liouville and Appell in ^ or can be derived from there (those 
belonging to Classes C and D). 

Another collection of Abel ODEs is found in the book by Murphy |^. After selecting those 
examples not having a constant invariant and for which a solution is shown in the book, we arrived 
at a set of nine ODEs, numbered in the book as: 78, 79, 80, 86, 275, 304, 345, 383 and 593. None of 
these ODEs represent an additional integrable class; their distribution among the classes discussed 
in this work is as follows 



Equations K.1.47, K.1.48 and K. 1.237 belong to Class C for infinitely many - however particular - values of one 
of the two parameters (see Eq.(E^); we don't know their solution for other values. 



Class 2 


Class 3 


Class B 


Class C 


Class D 


78, 80 


275 


86 


304, 383, 593 


79, 345 



Table 3. Classification for the non-constant invariant solvable Abel ODEs in Murphy's book. 



A wider collection of Abel ODEs than the one shown in Kamke's book is found in the book by 
Polyanin and Zaitsev ||T^. This book is rather new (1995) and covers a vast number of integrable 
ODE problems which we have not found in other books, hence making the examples attractive. 
On the other hand the Abel ODEs shown there are classified not according to their invariants but 
according to their form, and the origin of their solutions is not given. Apart from a main section 
consisting of four tables (82 Abel ODEs - all derived from four basic ones), the book contains other 
sections illustrating mappings between Abel and higher order ODEs. The quantity of examples is 
large and the computational routines we prepared for the equivalence problem are not yet covering 
in full the case in which the parameters of the class may assume symbolic values. As a result we still 
don't have a way to solve the equivalence problem for the whole set of integrable classes presented 
in |11|. Our analysis of these Abel ODEs of [0] is then still incomplete; consequently we restricted 
the presentation here to just a sample, constituted by the ODEs of the first of these four tables. 
These are 20 ODEs obtained from 



yy 



y = sx + Ax^' 



(45) 



by giving particular values to the parameters m and s {A is kept arbitrary). These ODEs appear 
in section 1.3.1 of |ll| under the numbers: 1, 2, 10, 16, 19, 22, 23, 26, 27, 30, 32, 33, 45, 46, 47, 
48, 53, 54, 55 and 56. We were not able to classify those numbered 27, 20, 48, 55 and 56. The 
distribution of the remaining ODEs, in the classes discussed in this work, is as follows: 



Constant invariant 


Class 1 


Class 2 


Class 3 


Class C 


Class D 


1, 2, 26 


23 


32 


33 


10, 19, 22, 45, 46, 47, 53, 54 


16 



Table 4. Classification for 15 of the 20 Abel ODE examples of Table 1.1 of juj 



5 New integrable Abel ODE classes derived from previous works 

Class "C" depending on one arbitrary parameter 

The form of the integrating factor studied by Abel actually leads to other integrable cases not 
mentioned in the original work Q . One of them is obtained by taking 6 = a in Eq. ( |3l| ) , resulting 
in the ODE family[3 



yy' - q'y 



(n+l) 







(46) 



where n / —1. The function q{x) and the parameter Ci can be removed as done with Eq. (|3^) , 
leading to 



y' = n[x- y^ - {n+ 1) 

which is turned exact by means of the integrating factor 



(47) 



(l + ((x2-x2'^)y-2x) y) 



n + l 
2 n 



(48) 



y " 



in Eq.@ is related to a in Eq.@ by n = -l/(2a + 1) 



1 



A simpler representative of this class is obtained by changing variables {y = u{t)t"-^ , x = t^-^}, 
then introducing a new parameter by means of n = 2;^, arriving at 

^ 2x ^ 2x ^ ' 

Taking into account Eq.(^8|), an implicit solution for this class is given by 

Class "D" depending on one arbitrary parameter 

In Q , Appell showed a series of integrable cases derived from the solutions to 



u' = A{u)+B{u)t (51) 



By changing variables {t = - — 4pT, u = x}, this ODE is transformed into the Abel ODE 



y' = -€-f4Vy' (52) 



B \B, 

where A and B are now functions of x. Any particular {A, i?} leading to a solvable case in Eq.([5l|) 
will then also lead to an integrable Abel ODE Eg. (|52|) . Among the choices for {A, 5} considered 
in - such that Eq.(^) results linear, homogeneous, or of Riccati type - only this mapping into 
Riccati type leads to something new. This case is obtained by taking 

A = ax^ + hx + c, B = ax^ + + 7 (53) 
The related Abel ODE family, depending on six parameters {a, 6, c, a, /3, 7}, is given by 

, 2 d ( ax'^ + bx + c \ 

ax^ + pa; + 7 ax \a x"^ + p x + j I 

and its solution could be expressed in terms of the solution to the Riccati ODE 

y' = {a + ax) y'^ + {b + P x) y + c + jx (55) 

However, we were not able to solve this Riccati ODE for arbitrary values of the six parameters 
involved and in |^ there is no indication of how that could be done. The alternative we then 
investigated is to consider the second order ODE obtained by replacing y = y' in Eq.(|55|). That 
ODE has two point symmetries if and only if a = 0. With these symmetries we were able to 
solve that second order ODE, and hence Eq.(^) when a = 0. Concerning the related Abel family 
Eq.(|5^ - now depending on five parameters - an appropriate change of variables of the form Eq.(§) 

|x = ^-^, y = v^n(t)| (56) 

followed by the introduction of a new parameter C by means of 

_ {(3^c + aj'^)a-ap-fb 
^ - p2 (57) 



transforms Eg. ([5^) into a simpler representative for the class 



l_iC^ (58) 

X x^ 

also showing that this class depends not on five but on one parameter^. It appeared of value to 
us also to determine the number of parameters on which Eq.(54) depends in the general case, that 
is before taking a = 0. For that purpose we searched for the appropriate changes of variables of 
the form Eq.(|^ which would remove as many as possible of these parameters, requiring that both 
the change and its inverse are finite. We then considered the branches which become infinite for 
some particular values of the parameters {a, 6, c, a, /5, 7} entering the transformations found. The 
results are summarized as follows. If all these parameters are different from zero, introducing new 
parameters {A, C, G} by means of 



a 



47 

A'^BC + l6-f^aA^B + pC + Ap^-f^ aA^B 
A'^B{/3^ + 4A^f 

C 



(59) 



C(2tA^ — B\ 

followed by changing variables {x = (^2_)_4^4)a^3g^ i U = u{t)A] in the six-parameter Eq.(|5J), one 
arrives at a 2-parameter representative for the same class 

G(Bx + x'^ -l)y^ 
y = — 7^ — + — ^ o-^^ (60) 

Now the case a = was already shown to lead to Eq. (|58|) , and all the other possible branches 
obtained from Eq.(^) by taking some of the other parameters equal to zero lead either to constant 
invariant families, or to members of the classes already discussed in this workp^ 
Three new classes not depending on parameters 

While analyzing the works ^, |^ and Kamke's examples, a large number of symbolic ex- 
periments were performed, sometimes leading to intermediate results which with a bit more work 
appeared to be new integrable classes by themselves. This happened three times, resulting in classes 
5, 6 and 7, for which representatives and solutions are given as follows: 
Class 5 

(2x + 3)(x + l)j/3 (5x + 8)y2 
y= 2^^ + 2.3 (61) 



■^^We note that in this process we have made two implicit assumptions: a 7^ and /9 7^ 0. To assure that the cases 



in Eq.(|54[) are covered by Eq.(|5q) we then also considered a = and P — separately, arriving at ODEs respectively 
members of the classes represented by Eq.(|5^) and Eq.(|35|). 

^''There is a special case, when & = 4 ^^"^'^4 , where the resulting Abel ODE can only be obtained from Eq.(|6ot) 
by taking appropriate limits. 



Solution: 



C, + , + iz- + l) dz = (62) 



4^ + 1 



where A = ^-^--^-A. 
Class 6 



Solution: 



^ (1 — X — x^) 
(x — 1)^ x^ (x — 1) 



i/' = -^r7^+ ' .V (63) 




xy (x — 1) / X — 1 + y 

—^dt is the exponential integral 



^(x-l)ye^ =0 (64) 



Solution: 



where A 



, (4x^ + 5x^ + 1)^3 (l-4x2)y 
^= +^ + 2x(x^ + l) 



Ci + 2 ,, +/ (z^ + l) ''^dz = Q (66) 

4 Q„,^2 



X — 2 yx'' — 3 yx'^ — y 



X (x + yx^ + y) 

6 Computer algebra routines, tests and performance 

The two itemized algorithms described in sections and ^ for solving the equivalence problem 
between two given Abel ODEs were implemented in Maple R5, in the framework of the ODEtools 



package [12|. The implementation consists of various routines, mainly accomplishing the following: 



1. determine whether a given Abel ODE belongs to one of the solvable classes described in the 
previous sections; in doing that, determine also the function F{t) entering Eq.(^) and the 
value of the parameters in the case of a parameterized class; 

2. use that information to determine the functions P{t) and Q{t) entering Eq.(^) and return a 
solution to the given ODE by means of changing variables in the solution available for the 
representative of the class (see the Appendix) . 

6.1 Representatives with simpler invariants for Classes A, C and D 

While preparing the computational routines being presented, we noticed that: if on the one hand 
the solutions to the representatives for the parameterized classes A, C and D (Eqs.(|33|, 58)) 



can be expressed in a relatively simple manner (see the Appendix), on the other hand, for each 
of these representatives, the form of the three invariants entering Eq.(p2[) is much simpler if an 
appropriate redefinition of the class parameter followed by a change of variables of the form Eq. (^) 
is performed. In turn, the complexity of these invariants is a relevant issue for a computer algebra 
implementation since simpler invariants lead to simpler GCD and resultant computations. 



In the case of Class A (Eq.(|33|)), redefining the parameter as a = 5/36 — k/3 - where k is the new 
parameter - and changing variables {x = ^2t-~i^ ' ^ ~ ^ Vi2^<-6 ^ ^^^"^ ^° class representative 

, ^ {{3k-2x-x^) y^-9y^-9y) 
^ 9(2x-l) ^ ' 



and the first invariant S5/S3 for it is 

(9k^ + (SOx^ - 6x + 3) K- + 5x'^ + 7 x'^f 
9 (k + x2)^ 

as opposed to 



(68) 



(27 g^x^ - (72 x'^ + 270 + 15 x^) a + 54 + 36 + 135 + 15x^ + 2 x^)^ 
x4(2x4-9ax4 + 9x2 + 27)^ 

which is the same first invariant but calculated for Eq.(|33|). Actually a measure of the Maple 
computational length of the three invariants [sf/sl, s^sj/s^, S5S7/S3] entering Eg. ([2^) shows the 
values [99, 215, 223] when calculated on Eq.(67) above and [146, 354, 356] when calculated on 



Eq.(33). Also, the degrees in x of both numerator and denominator of Eq.(|68D are lower than the 
corresponding degrees of Eq. (|69|) . 

A similar situation happens with Class C (Eq.(|49|)), where redefining the parameter as a = 



and changing variables | x = f") (■^ *+^) i ^ y _ t u{t)K{{3t+4)K+3t+4) \ j^^^^ ^ class represen- 

^ (K-rL)i ^(k+1)(3— ft)(3 1+4) t J 

tative a bit more complicated than Eq.(|49|): 

' A ff , Q^ 3 , r 2 {6KX + 5K + 3)y 

y =4 ((k-2 x + K-3 Kxy +6y — (70) 

(3x + 4) Kx 

for which, however, the computational length of the three invariants [sf/sl, s^sj/s^^ sssy/sl] is 
[110, 273, 305] as opposed to [173, 382, 412] when calculated for Eq.(|4^). The degrees with respect 
to X of the numerators and denominators of the invariants of Eq.([70|) are also lower than those of 
the invariants of Eq. (|49|) . 

The same reduction in the complexity of the invariants is achieved for class D (Eq.(p8|)) by redefining 



the parameter via C = | \/k and changing variables {x = J^-^ ^ V = y^2 (1 — t^/t) y/Ku{t)}, 

leading to 

I {2x^K-2)Ky^ -3y'^K + Kxy 

y = 



for which the first invariant, S5/S3, is given byp^ 

((2x^ + 15x3) ^2 _ (43.2 _^ 153. _ 9) «;4.2)' 



4k^ (kx^ — X + 1)^ 

as opposed to 



(72) 



(C + x^)^ (2 + 8 x'^C^ + (12 x^ + 15 x2) C2 + 8 x^C + 2 x^ - 15 x^ + 9 x^)^ 
729 z (73 1 

{2C^ + 6C^x'^ + {9x^ + 6x^)C + 2x^ -dx'^f 



The case k = is treated separately. 



which is the same first invariant but calculated for the class representative Eg. (pq) . In Eq.(72), not 
only the degrees with respect to x but also those with respect to the class parameter k are lower 
than the corresponding degrees in Eq.(j7^). 

Moreover, for class D this change in the representative of the class also fixes a problem: the 
representative Eq.(58) is itself invariant under C — > —C followed by the change of variables 



iC_ 



y = iu{t) 



(74) 



where i is the imaginary unit, and hence Eq.(73) is invariant under {C — > —C, x — > iC/x}. 
Consequently, if {C, F{t), P{t), Q{t)} is a solution to the equivalence problem between a given 
ODE and Eq.(^8|), then {— C, -p^, iPit), Q{t)} is also a solution; this fact invalidates the use of 
the interpolation method described in sec. with Eq.(|5^) since that method can only be used 
when the interpolated solution C(a) is unique. 



6.2 Installation 

The programs being presented have been written as one more step in the development of the 
ODEtools Maple package [12| and hence are integrated to it and not distributed separately. To 
install the new Abel routines then what is necessary is to install ODEtools to run in the Maple 
environment by putting the related libraries (two files - maple. ind and maple. lib) in any directory 
- say DEJibraries-directory - then opening Maple, and adding that directory to Maple's libname 
variable via 



> libname := DE_libraries_directory , libname: 



where '>' is the Maple prompt. This instruction automatically updates Maple's dsolve subroutines 
to make use of the new routines for solving Abel ODEs - no further steps are required. In this 
way, the new Abel ODE routines are automatically used by Maple's dsolve when the input ODE 
is of Abel type, as well as to solve higher order ODEs when they can be reduced to first order Abel 
ODEs members of the classes discussed in sec. ^ and ^ 

Apart from this integration with dsolve, it is also possible to try just the routines being pre- 
sented by giving to dsolve the extra argument [Abel] . For example (Kamke's ODE 37) 

> ode[37] := dif f (y (x) ,x) -y (x) "3-a*exp(x) *y (x) "2 = 0; 

odesT ■=y' -y^ - ae'y'^ = 



> dsolve (ode [37] , [Abel]); 

These implicit answers can also be tested in the Maple worksheet, interactively, using the standard 
Maple o detest command. 

Due to the intrinsic complexity of Abel equations, the solution for most of the solvable classes is 
expressed in implicit form and in terms of elliptic integrals and special or hyper geometric functions. 
Then, to save the time Maple spends in trying to "integrate" these integrals or to "invert" these 
algebraic expressions it is frequently convenient to call dsolve with the optional extra arguments 



'useint, implicit', meaning: use Maple's inert Int and return the solution directly in implicit 
formp^. For example, for Kamke's 185 



> ode[185] := x"7*dif f (y (x) ,x)+2* (x"2+l) *y (x) "3+5*x"3*y (x) "2 = 

odei85 := x^y' + 2 (x^ + l) + 5 = 
the solution obtained using these extra arguments is: 

> dsolve (ode [185] , [Abel], useint, implicit); 

C, + ^=J== + l/2 j\z- + l)-"'dz = ^ (76) 

where A = + The explicit computation of the integral above leads to a complicated 
expression with hyper geometric functions somehow obscuring the structure of the solution. 

Concerning Maple's difficulty in solving the problem when the class parameter C involves radicals 



(see end of sec. 3.1), we implemented an environment variable controlling how hard the routines 
will work, in order to avoid exhausting system resources unless a hard trial is specifically requested. 
This environment variable is _Env_odsolve_Abel_try_hard, it can be assigned a positive integer 
from 1 to 5, and by default it is assigned to 4; the meaning of the possible values is as follows: 

• if _Env_odsolve_Abel_try_hard = 1 then the algorithm for parameterized Abel classes is 
disabled, and for non-parameterized classes only a restricted equivalence using {x = t,y = 
P{t) u + Q}, that is: with F{t) = t in Eq.(|), is tried; 

• if _Env_odsolve_Abel_try_hard = 2 then for non-parameterized classes a full equivalence 
using Eq.(§) is tried; 

• if _Env_odsolve_Abel_try_hard = 3 then the algorithm for parameterized Abel classes is 
enabled but only for numerical solutions for the class parameter C; 

• if _Env_odsolve_Abel_try_hard = 4 then for parameterized Abel classes both a numerical 
or a rational interpolation solution for the class parameter C are tried (symbolic variables in 
the coefficients are allowed but solutions with radicals are not computed) ; 



if _Env_odsolve_Abel_try_hard = 5 then the algorithm for parameterized Abel classes will 
also compute solutions involving radicals for the class parameter C. 



Finally, these routines for Abel ODEs were programmed to provide extensive run-time infor- 
mation on the computations being performed through the Maple standard userinfo & infolevel 
scheme. For example, Kamke's ODE 43 belongs to class B (see Table 2.) and both the way how 
this is determined and the explicit values found for F, P, Q, and the parameter entering Eq.(|37|) 
can be seen by setting the infolevel as followsp^: 

> infolevel [dsolve] := 4; 

> ode [43] : = 



^^The option implicit is standard in Maple's dsolve and the implementation of the option useint comes with 
the routines for Abel ODEs presented here. 

We kept few lines to illustrate the userinfo feature, and represented the missing ones by 



ode43 := y' + (Sax^ + 4a^a; + 6j + = 
> dsolveCode [43] , [Abel], implicit); 

The relative invariant s3 is: -3*a*x~2+b-2*x~3 

The first absolute invariant s5"3/s3"5 is: 108*(12*a"2*x"3-8*a*x*b + .. 

The second absolute invariant s3*s7/s5''2 is: 1/3* (3*a*x''2-b+2*x''3) * .. 
The third absolute invariant s5*s7/s3''4 is: 9* (372*a"3*x"4+450*a"2*x"5 

-> ====================================== 

-> ...checking Abel class B (by Liouville) 

Trying a = 
Trying a = 1 
Trying b = 
Trying b = 1 

-> Step 1: checking for a disqualifying factor on F after evaluating x 
Trying x = 

*** No disqualifying factor on F was found *** 

-> Step 2: calculating resultants to eliminate F and get candidates for 

*** Candidates for C are [1, 4, 1/4], *** 

-> Step 3: looking for a solution F depending on x 



C = 1/4 leads to the solutions [{F = -l-3/2*x}] 
Interpolated candidate for the class parameter C is: C = 1/4 
General testing of the candidate C = 1/4 with arbitrary b 
Interpolation is still incomplete; trying next value of b 



Trying b = 2 

General testing of the candidate C = -3/4*b+l with arbitrary b 

C = -3/4*b+l leads to the solutions [{F = -l-3/2*x}] 
General test of C = -3/4*b+l passed DK; 
interpolation for b in this level is complete. 

General testing of the candidate C = (-3/4*b+a"3)/a"3 with arbitrary a 

C = (-3/4*b+a~3)/a~3 leads to the solutions [{F = -l/2*(3*x+2*a)/a}] 
General test of C = (-3/4*b+a"3)/a~3 passed OK; 
interpolation for b in this level is complete. 



Value of the Class parameter solving the problem is : 

C = l/4*(-3*b+4*a-3)/a-3 
Inverse of the transformation solving the problem is: 



{t = -l/2*(3*x+2*a)/a, u(t) 



= -2/3*a"2*y(x)} 



Solution: 

Cl + ^TT^ -^ / = (77) 

((^^ - ^1) KAi, -VM) + VMKAi + 1, -VM)) 

where Ai = ^ ^^^^5^' ^2 = + ^^^J^ " 2^' and I and K are respectively the modified 

Bessel functions of first and second kind. 



6.3 Tests and Performance 

The idea was to test these computational routines to confirm the correctness of the returned so- 
lutions as well as to obtain the classification presented in the previous sections for solvable Abel 
ODEs. The first testing arena was the 69 Abel examples found in Kamke plus the 9 solvable exam- 
ples with non-constant invariant from Murphy's book, plus the first 20 Abel ODE examples from 
[0] mentioned at the end of sec. |^ - totaling 98 Abel ODE examples. The routines passed these 
tests - the solutions obtained were confirmed to be correct using other symbolic computation tools 
interactively - and the resulting classification is that shown in Tables 1, 2, 3 and 4 of sec. 

The aforementioned test however involves only 40 "non-constant invariant and solvable" Abel 
ODE examples, and does not fully test the new routines. We have then set up a more thorough 
test, which can be taken as a test-suite for other computer algebra implementations of methods for 
solving Abel ODEs. The ideas behind this additional test-suite are summarized as follows: 

• Take the representatives of the seven Abel class not depending on parameters, discussed in 
sec. Q and |5|, and generate with each four more Abel ODEs of the same class by applying the 
different types of transformations: 



1. The general case with F, P and Q arbitrary: 

tn := {x = y{x) = P{t)u{t) + Q{t)} (78) 

2. Rational transformation with symbolic coefficients a and b: 

tra :={x = ^ + y{x) = ^ + l} (79) 



3. Non-rational transformation involving an elementary function and symbols: 

trs := {x = e* + 1 + -, y{x) = u{t) + 1} (80) 



4. Non-rational transformation involving abstract powers and symbols: 

tr^:={x = ^ + i, y{x) = uit) + l} (81) 



For example, the first of these four Abel ODEs generated from Eq.(p^) is given by: 



F'p2(_3 + i?)y3 ((3F-9)Q- 10)PF'y2 



8F 



+ 



(82) 



+ 



■((3F-9) Q2_20Q-3) F' 



8F 




and the other three are obtained from this one by replacing F, P and Q by the values implied by 
Eqs. ([79|), (^) and (|8ll). The scheme just outlined generates 28 more solvable Abel ODE examples 
with non-constant invariant (4 per class), and suffices for testing the solving of classes without 
parameters. 

Concerning classes with parameters: 

• Take the representatives of the four Abel classes depending on parameters discussed in sec. 
^ and 6.1, generate with each one four more Abel ODEs of the same class by applying the trans- 
formations Eg. ([raj) and Eq.(^), preceded by replacing the single parameter C entering each class 
representative by C = 2 and C = a/ {3. 

This increases by 4 x 2 x 2 = 16 more ODE examples, and the solving of these examples tests 
both the scheme for numeric values of the parameter C as well as the case where C is a rational 
function of other symbols entering the given ODE. 

The time spent by the routines being presented in solving all these 28 -|- 16 = 44 additional 
Abel ODEs with non-constant invariant, is summarized in Tables 5, 6 and 7 as follows: 



Clas 






transf. (i) 


transf. (ii) 


transf. (iii) 


transf. (iv) 


1; Eq.( 


34 


) 


0.255 sec. 


0.879 sec. 


0.687 sec. 


1.256 sec. 


2; Eq.( 


35 

i 


) 


1.530 sec. 


1.900 sec. 


1.998 sec. 


2.620 sec. 


3; Eq.( 


) 


0.277 sec. 


6.760 sec. 


5.623 sec. 


11.635 sec. 


4; Eq.( 


) 


0.456 sec. 


7.126 sec. 


3.941 sec. 


15.511 sec. 


5; Eq.( 


3] 

1 

35 


) 


0.610 sec. 


2.505 sec. 


3.082 sec. 


7.877 sec. 


6; Eq.( 


) 


1.379 sec. 


120.952 sec. 


45.565 sec. 


282.103 sec. 


7; Eq.( 


) 


0.896 sec. 


15.441 sec. 


24.586 sec. 


183.557 sec. 



Table 5. Timings for 28 Abel ODEs with non-constant invariant - 7 classes free of parameters. 



Class parameter C = 2 


Class 


transf. (ii) 


transf. (iv) 


A; Eq.( 


37 
37 


) 


25.546 sec. 


41.304 sec. 


B; Eq.( 


) 


29.396 sec. 


54.171 sec. 


C; Eq.( 


70 


) 


18.648 sec. 


41.425 sec. 


D; Eq.( 


71 


) 


36.641 sec. 


75.962 sec. 


Table 6. Timings for 8 Abel ODEs with non-constant invariant - "numeric" class parameter. 


Class parameter C = ^ 


Class 


transf. (ii) 


transf. (iv) 


A; Eq.f 


37 
37 


) 


67.076 sec. 


156.425 sec. 


B; Eq.( 


} 


174.751 sec. 


491.210 sec. 


C; Eq.( 


70 


) 


73.940 sec. 


174.761 sec. 


D; Eq.( 


71 


) 


88.365 sec. 


177.314 sec. 



Table 7. Timings for 8 Abel ODEs with non-constant invariant - "symbolic" class parameter. 



Concerning a comparison of performances between the new routines and those available in other 
computer algebra systems (CAS), this appeared to us not justified in this case: roughly speaking 



none of these CAS return solutions for Abel ODEs with non-constant invariant. More precisely, 
from Mathematica 3.0, Macsyma 2.7, Maple 5.1 and Mupad 4.0, all of them failedp^ in solving any 
of the 18 Kamke's examples shown in Table 2. (and hence in solving any of the 44 ODEs of Tables 
5, 6 and 7), except for Kamke's example 235 - it is an inverse linear ODE - and anyway none of 
them solved it after transforming it to first kind format. 

6.4 Performance with the 1*^* order ODE examples from Kamke 

Although the main purpose of this paper is to present a computational scheme for finding solutions 
to Abel ODEs, it is interesting to see how odsolve - the ODE-solver of the ODEtools Maple package 
[ p!^ ] - performs with the addition of these new routines. The performance with all of Kamke's 553 
solvable examples^ after incorporating the computational routines presented in this paper is: 97% 
are solved. This performance is summarized as follows 



Degree in y' 


ODEs 


Solved 


Average time 


solved 


fail 


1 


350 


337 


3.2 sec. 


12.9 sec. 


2 


145 


140 


8.8 sec. 


61.1 sec. 


3 


27 


26 


7.2 sec. 


17 sec. 


higher 


31 


30 


13.4 sec. 


25.2 sec. 


Totah 


553 


533 


~ 6 sec. 


^ 20 sec. 



Table 8. Kamke's first order ODEs, solved by odsolve; 97% 



The number and classification of Kamke's 1*^* order ODEs still not solved by odsolve is now: 



Class 


Kamke's numbering 


rational 


452 


Riccati 


25 


NONE 


80, 81, 83, 87, 121, 128, 340, 367, 395, 460, 506, 510, 543, 572 



Table 9. Kamke's f order solvable ODEs for which odsolve fails: 3% 



where the Abel ODEs numbered in Kamke's book as 47, 48, 205, 206, 237 and 265 not presented 
in the tables above are known to be solvable only for specific values of their parameters - and for 
these values odsolve succeeds - and the ODEs 234 and 253 were not included since their solutions 
are not shown in the book or known to us. 

7 Conclusions 

In this paper, a first classification, according to invariant theory, of solvable non-constant invariant 
Abel ODEs found in the literature, was presented. Also, a set of new solvable classes, depending 
on one or no parameters, derived from the analysis of the works by Abel, Liouville and Appell 

^, ^, was shown. Computer algebra routines were then developed, in the framework of the 
Maple ODEtools package, to solve - in principle - any member of these classes by solving its related 
equivalence problem. The result is a concrete new tool for solving Abel type ODEs fully integrated 
with Maple's ODE-solver dsolve. 

The classification shown has had the intention of giving a first step towards organizing in a 
single reference the integrable cases widely scattered throughout the literature. The derivation of 

■^^For Kamke's ODE 257, Macsyma (2.7) returns a wrong answer in terms of y' . 

^**We classified as unsolvable in general Kamke's examples 50, 55, 56, 74, 79, 82, 202, 219, 250, 269, 331, 370, 461, 
503 and 576. 



new solvable parameterized classes from the works by Abel and Appell in the 19*'^ century (Classes 
"C" and "D") also showed that valuable information can still be obtained from these old papers. In 
fact, from Tables 1, 2, 3 and 4 in sec. ^ the larger number of integrable cases found in textbooks are 
particular members of this Class "C" (Eq.(^9|)) - an integrable class derived from a case somehow 
disregarded in Abel's Memoires [|[. 

As for the computer routines, the implementation presented here for solving the equivalence 
problem for parameterized classes proved to be a valuable tool in most of the Abel ODE examples 
we were able to collect. The Abel ODE routines here presented were recently integrated to the 
Maple computer algebra system and are already part of the upcoming Maple release 6. 

On the other hand, we note the intrinsic limitation of this Abel ODE problem: most of the 
solutions can only be obtained in implicit form and in terms of quadratures; in turn, these integrals 
are usually elliptic integrals so that they cannot be expressed using elementary functions. Also the 
classification of integrable cases presented here is incomplete in that it is missing - at least - a more 
thorough analysis of the integrable cases presented in |11|. We are presently working on this topic 
and expect to succeed in obtaining reportable results in the near future. 
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Appendix A|^ 



Class 


Representative equation and solution 


1 


1 3y^~3y-x ^ , ^ ('^ ^ +(^V -3by+27)x+4y -iy ) 
y x(8y-9) ' ^1 1 {x^+2x{y^~3y)+y^f 


2 


xAilx^-^]+Alll,x'^-^] 


3 




4 


y' = y'-^J/', Ci + iet-^-Ei(l,x-i)=0 


5 


/_ {2x+3)(x + l)y' (5x+8)y^ 
y 2x'^ ' 2x-^ ' 

4/ (.+ 1)2 ^ ^ ^ 4 


6 


- vl^ + , Ci Ei f 1, ''Y'-^) + (^-^)y^ f-'> - 

a;2(j._i-)2 x-'{x — l)'' ' ' a:a(a: — 1) y x — l + y 


7 


y 2x^ 1 y 1 2a:(:r2 + l) 

•' ^ ' x{x + yx^+y) 


A 


y' = {ax+l + ^)y' + y'', 

Ci + ^x e-P (j^l^^^) - /^exp (Z,,^^) = 


B 


y' = 2 {x'~a')y' + 2{x + l)y', 

{a + x)K(^a,-^x^ + ij-a^'j + ^x2 + -^-a^ K + a , - ^a:^ + i -q2 j 


C 


j/' = ^^^^ + («-i)y^-ff- c. + f (i-(i-x)^)''"-2/'^(i-.^)^d. = o 


D 


, _ y^ (a + x'^)y^ 

^ (.+l)M(-f-|.i.i(.-f-l)^)+(|-.^)M(-f + i.i,i(.-f-l)^) ^ 



Representative ODEs and their solutions for the Abel ODE classes presented in this work. 



The solution shown for the representative of class D is not valid when a is an integer, or when 2a is a positive 
integer. In those cases, the solution of the associated Riccati equation Eq.(^^ takes many different forms depending 

on the value of a, which we found inconvenient to present here. Ei(n,3;) — / dt is the exponential integral, 

Ai(a::) and Bi(a;) are the Airy wave functions, K(a;) and \{x) are the modified Bessel functions of the first and second 
kinds, respectively, and M(a;) and W(x) are the Whittaker functions. 



